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SUMMARY 

A formal procedure for the probabilistic design assessment of a composite structure is described. The 
uncertainties in all aspects of a composite structure (constituent material properties, fabrication variables, struc- 
tural geometry, and service environments, etc.), which result in the uncertain behavior in the composite struc- 
tural responses, are included in the assessment. The probabilistic assessment consists of (1) design criteria, 

(2) modeling of composite structures and uncertainties, (3) simulation methods, and (4) the decision-making 
process. A sample case is presented to illustrate the formal procedure and to demonstrate that composite struc- 
tural designs can be probabilistically assessed with accuracy and efficiency. 


INTRODUCTION 

Composite materials are widely used in modem structures for high performance and reliability. However, 
because these structures usually operate in hostile and random service environments, it is difficult to predict 
the structural performance. In addition, experiments show that the composite structural behavior exhibits wide 
scatter as a result of the inherent uncertainties in design variables. The design variables, known as primitive 
variables, include the fiber and matrix material properties at the constituent level; fiber and void volume ratios, 
ply misalignment and ply thickness for the fabrication process; and random structure size, loadings, and 
temperature. 

The scatter in the structural behavior cannot be computationally simulated by the traditional deterministic 
methods that use a safety factor to account for uncertain structural behavior; thus, the structural reliability can- 
not be discerned. A probabilistic design methodology is needed to accurately determine the structural reliability 
of a composite structure. In the past, Monte Carlo methods have been widely used for probabilistic composite 
structural analysis (ref. 1). However, these methods are computationally intensive and can be used for verifi- 
cation purposes only. Likewise, when perturbation techniques (ref. 2) are used to account for uncertainties, only 
the first few statistical moments of the structural responses — and not the cumulative distribution function 
(CDF) — are obtained. 

NASA Lewis Research Center has developed a formal methodology to efficiently and accurately quantify 
the scatter in the composite structural response and to assess the composite structural design, while accounting 
for the uncertainties at all composite levels (constituent, ply, laminate, and structure) (fig.l). This methodology, 
which integrates micro- and macrocomposite mechanics and laminate theories, finite element methods, and prob- 
ability algorithms, was implemented through the computer code IP ACS (Integrated Probabilistic Assessment of 
Composite Structures) (fig. 2) (ref. 3). EPACS is used to assess composite structures probabilistically for all 
types of structural performances such as instability, clearance, damage initiation, delamination, microbuckling, 


fiber crushing, and resonance damage. Since IPACS uses a special probability algorithm FPI (fast probability 
integrator) (ref. 4) instead of the conventional Monte Carlo simulation, tremendous computational time can be 
saved (ref. 5). Therefore, a probabilistic composite structural analysis, which cannot be done traditionally, 
becomes desirable especially for large structures with many uncertain variables. A typical case is analyzed 
herein to demonstrate the code IPACS for the probabilistic assessment of composite structures and to illustrate 
the formal design assessment methodology. 


FUNDAMENTAL APPROACH 

The fundamental approach for the probabilistic assessment is as follows: 

(1) Identify all possible uncertain variables at all composite scale levels. 

(2) Determine the probabilistic distribution function (PDF) for each variable. 

( 3 ) Process all random variables through an analyzer that consists of micro- and macrocomposite mechanics 
and laminate theories, structural mechanics, and probability theories. 

(4) Extract useful information from the output of the analyzer and check against the probabilistic design 
criteria. 

The uncertainties in a composite structural design can originate at different composite scale levels. At the 
constituent level, the material properties for the fiber and matrix are the major sources of uncertainties. The 
primitive variables are defined in the appendix and their typical values are listed in table I. At all stages of the 
fabrication process, the fabrication variables such as fiber volume ratio, void volume ratio, ply misalignment, and 
ply thickness show considerable scatter. At the structure level, variation of the geometry during the assembly stage, 
uncertain boundary conditions, and random thermal-mechanical loads contribute significantly to the scatter in the 
composite structural response. 

Once the uncertainties in the primitive variables are identified, micro- and macrocomposite mechanics and 
laminate theories are used to propagate them in the constituent material properties and in the fabrication vari- 
ables from the lower composite levels (ply) to the higher composite levels (laminate). The scatter in the struc- 
tural response is then computationally simulated through structural mechanics (or finite element methods) and 
through probability theories (Monte Carlo simulation or FPI) to account for uncertainties in the laminate mate- 
rial properties, geometry, boundary conditions, and thermal-mechanical loads. 

Probabilistic design assessment includes probabilistic design criteria that are based on considerations such 
as safety, performance, economics, and other requirements that could introduce uncertainties in a specific design. 
These factors contribute to the uncertainties of the structural response and provide inclusive information for 
judicious decisions on design acceptance or rejection. 


PROBABILISTIC COMPOSITE STRUCTURAL DESIGN ASSESSMENT 

This report describes a probabilistic design assessment methodology for composite structures with an 
acceptable or preassigned reliability. The details of the assessment are described in the following three steps. 
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Step 1 : Probabilistic Design Criteria Setup 

A typical design criterion can be stated as follows: The probability of a failure event should be less than an 
acceptable value, say 10' 3 . A failure event occurs when a structural response is less than the allowable response. 
This probability is defined as the failure probability. The allowable response divides the possible response 
domain into safe and failure regions as shown in figure 3. The predicted failure probability is the area under the 
probability density function in the failure region. The critical response (fig. 3) is determined by IPACS such that 
the probability of a response exceeding this critical value is in the safe region. When the critical response falls 
within the safe region, the design is acceptable. When the critical response falls within the failure region, the 
design is unacceptable and requires a redesign. Sample probabilistic design criteria for the various failure modes 
are described as follows: 

(1) Instability— The probability that the buckling load is smaller than the design load should be less than 

10 ' 3 - 


(2) Clearance — The probability that the nodal displacement is greater than the allowable tolerance should be 
less than 10' 3 . 

(3) Resonance avoidance— The probability that the natural frequency is greater that its upper bound should 
be less than 10' 3 ' 

(4) Delamination — The probability of delamination occurrence should be less than 10' ‘ 


Step 2: Probabilistic Simulation Using IPACS 

IPACS integrates several NASA in-house computer codes developed in recent years such as COBSTRAN 
(Composite Blade Structural Analyzer) (ref.6), PICAN (Probabilistic Integrated Composite Analyzer) (ref. 7) and 
MHOST (Marc Hot Section Technology) (ref. 8). COBSTRAN is a dedicated finite element model generator for 
structural analysis of composite structures. PICAN uses ICAN (Integrated Composite Analyzer) computer code 
(ref. 9) for composite mechanics. This code has evolved over the last 20 years and has been verified with 
experimental data for all aspects of composites. PICAN enables the computation of the perturbed and probabilis- 
tic composite material properties at the ply and laminate levels. MHOST performs structural analyses using veri- 
fied finite element methods. These analyses determine the perturbed and probabilistic structural response at 
global, laminate, and ply levels. PICAN and MHOST share the FPI module (ref. 4) for the application of the 
fast probability integration algorithm to obtain cumulative distribution functions (CDF's) of the material pro 
perries and the structural responses. 

FPI is an approximate technique for the probabilistic analysis of the structural performance and has the 
major advantage of speed. FPI techniques are orders of magnitude more efficient than Monte Carlo simularion 
methods. This is especially true in the tails regions of the distribution; that is, at very high or low probabilities 
since the FPI solution time is independent of the probability level. Conversely, in Monte Carlo simulation 
methods, the computational time increases with very high or low probability levels. Also, FPI allows evaluation 
of information that describes the relative importance of each random variable. These sensitivity factors can be a 
valuable aid in optimizing a design. 

In IPACS, the probabilistic assessment of composite structures starts with the identification of uncertain 
primitive variables at constituent and ply levels. These variables are then selectively perturbed several times to 
create a data base. The data base is used to establish the relationship between the desired structural response (or 
the desired material property) and the primitive variables. For every given perturbed variable, micromechanics is 
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applied to determine the corresponding perturbed mechanical properties at the ply and laminate levels. Laminate 
theory is then used to ascertain the resultant force-moment-strain-curvature relationship. With this relationship at 
the laminate level, a finite element perturbation analysis is performed to find the structural response that 
corresponds to the selectively perturbed primitive variables. This process is repeated until enough data are 
generated and the proper relationship between structural response and primitive variables can be established 
through a numerical procedure. 

Given the probabilistic distributions of primitive variables and a numerically determined relationship 
between them and the structural response, FPI is applied. For every discrete response value, a corresponding 
cumulative probability can be computed quickly by FPI. This process is repeated until the CDF can be appro 
priately represented. The probabilistic material properties at ply and laminate levels are also computed in this 
same way. The output information from FPI for a given structural response includes its discrete CDF values, the 
coefficients for the PDF that was used for the uncertainties in the primitive variables, and the variables’ 
sensitivity factors to the structural response. It is important to restate that the process consists of multiple 
deterministic analyses using verified structural and composite mechanics computer codes. Therefore, the 
probabilistic simulations are mechanistically accurate and reflect the uncertainties of the primitive variables in 
the uncertainties of the structural response. 


Step 3: Decision Making and Redesign 

IP ACS simulates the PDF of a given structural response, such as buckling load, displacement, local stress, 
local strength, vibration frequencies, and fatigue life. The probability of a design violation for each criterion can 
be calculated with these PDF functions. When the failure probability is greater than the acceptable value, say 
10' 3 , the composite structural design should be rejected. To redesign a composite structure, one can use the 
sensitivity factors from the IPACS analysis. Sensitivity factors rank the random variables based on their 
contribution to this failure probability. Therefore, a redesign will be guided by this information with manufactur- 
ing control of the mean and the standard deviation (stdv) of the appropriate random variables. 


DEMONSTRATION CASE AND DISCUSSION 

A stiffened composite cylindrical pipe is probabilistically assessed against probabilistic design criteria. The 
cylindrical pipe is 2 ft in diameter and 20 ft long (fig. 4). The structure is modeled by 588 four-noded shell 
elements and 600 active nodes (6 degrees of freedom per node). The composite pipe consists of the skin, three 
horizontal circumferential frames, and four vertical stringers. The laminate configurations for the skin, frames, 
and stringers are [±45/0^5/0 2 /±4S76m\, [0 24 ] and [0 24 ], respectively. The pipe is assumed to be supported at 
one end by a set of translational and torsional spring constants and free at the other end. When the spring con- 
stant approaches infinity, a completely fixed boundary condition is simulated. When the spring constants are set 
to zero, a free boundary condition is simulated. For a given set of spring constants, a partially fixed boundary 
condition is modeled. The pipe is subjected to axial (F x ) and lateral (Fy) loads as well as torsional moments 
(M^) at its free end (fig. 5). 

The uncertain variables are identified at the constituent, ply, and structure levels. At the constituent level, 

17 material properties for the graphite fiber and 12 for the epoxy matrix of the skin, frames, and stringers are 
assumed to be uncertain variables. Their probability distribution types and associated parameters are listed in 
table I. At the ply level, the fabrication variables (fiber and void volume ratio, ply orientation, and ply thick- 
ness) are treated as random variables. Their statistics are shown in table II. At the structure level, spring 
constants that simulate a partially fixed boundary condition are assigned by a probability distribution as are the 
loading conditions. Their statistics are shown in table III. 
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In the following paragraphs, the composite pipe is assessed or checked against two design criteria, clearance 
and delamination; and the results are discussed. 


Clearance Assessment 

The clearance criterion is violated when the displacement at the free end in the lateral direction is greater 
than the allowable value. In this assessment, acceptable failure probability is chosen to be 10 . From the static 
analysis, the probabilistic displacement at the free end in the lateral direction was simulated as shown in fig- 
ure 6. The critical displacement corresponding to 10‘ 3 failure probability is 5.2 in. If the allowable displacement 
is 6 in., then the critical displacement falls in the safe region, and the clearance criterion is satisfied. If the 
allowable displacement is 5 in., then the critical displacement falls in the failure region. The clearance design 
criterion is violated and the pipe needs to be redesigned. From the IP ACS sensitivity analysis, the fiber 
modulus, fiber volume ratio, ply thickness of the skin, and random loads in lateral direction have the most 
significant contribution to the failure probability (fig. 7). 


Delamination Assessment 


Delamination occurs when the ply stress is greater than the ply delammation strength. From the IP ACS 
analysis, the relationship between the ply stress S pl and the independent random variables X is numerically 
determined as shown in equation (1). 


N 

Spl = ao + l 


a, X 


N 

*1 

i=l 


x i2 


(i) 


where a^, a ; and b; are constants; N is number of independent random variables. The ply delamination strength 
S DL (ref. 9) is shown in equation (2). 


Sdl — 10 Sj + 2.5 S mT A 


( 2 ) 


S,= 




A 


(3) 


and 


A = 1 - 


4wr 

> n 1 - fvr] 


( 4 ) 


where fvr and wr are fiber volume ratio and void volume ratio; G^ and G m are fiber and matrix shear 
modulus; S mT and S mS are the matrix tensile and shear strength. A limit state function (LSF) is defined as 

LSF = S DL - S pl 
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If the LSF is less than 0, it indicates that the ply stress is greater than the ply delamination strength. There- 
fore, delamination will occur. The probability of LSF < 0, computed by FPI, is 0.0008, which is smaller than 
the acceptable failure probability (10' 3 ). From the sensitivity analysis, the eight most influential random vari- 
ables that contribute to the failure probability are identified in figure 8. The matrix shear strength is the most 
important (sensitivity factor about 0.8) followed by the ply thickness of the composite skin (sensitivity factor 
about 0.4). If the acceptable failure probability is reduced to 1CT 4 , redesign is necessary. The redesign can be 
achieved most efficiently by manufacturing tolerance control or by controlling the mean or scatter of the signifi- 
cant primitive variables. For example, if the coefficient of variation (scatter) of the ply shear strength is reduced 
from 5 to 4 percent, the failure probability is reduced to 0.0003. However, if the scatter of the ply thickness is 
also reduced from 5 to 4 percent, the failure probability can only be reduced to 0.0006. This demonstrates that 
the failure probability can be reduced more effectively by reducing the scatter of the ply shear strength, the most 
important random variable. 

Any other structural response can be similarly evaluated. These two examples demonstrate how the probabi- 
listic design assessment is performed by using IPACS. 


CONCLUDING REMARKS 

A formal methodology is described in this report for the probabilistic design assessment of composite 
structures. This methodology, integrating micro- and macrocomposite theory, structural mechanics (finite 
element methods), and probability algorithms, performs a probabilistic assessment of composite structural 
designs considering all identifiable uncertain variables at all composite levels. Composite structural designs can 
be assessed against specific probabilistic design criteria demonstrating that such designs can be computationally 
assessed by using the probabilistic computer code IPACS. Information for an efficient design can also be 
obtained. Specifically, for the demonstration case, the uncertainty range in the end displacement was between 1 
and 3 percent of the pipe length and was most sensitive to the uncertainties in the skin-related primitive 
variables. Conversely, probable delamination failure was most sensitive to the shear strength of the skin. 
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APPENDIX— SYMBOLS 


Cf 

C, 

D„ 


m 


E m 

E f22 

E m 

E x 

G fl2 

G f23 


K. 


cTR 


K cTO 

Kf„ 

K f22 

Kf33 

M xx 

N f 

S DL 

S fC 

SfT 

S mC 

S mS 

S mT 

S P' 

stdv 

*psk 

^pst 

wr 

X 


«fll 

a i22 

Pm 

% 

V fl2 

v f23 

V m 

Pf 

Pm 


fiber heat capacity, Btu/in. °F 

matrix heat capacity, Btu/in. °F 

matrix diffusivity, in. 3 /sec 

filament equivalent diameter, in. 

fiber modulus in longitudinal direction, Mpsi 

fiber modulus in transverse direction, Mpsi 

matrix elastic modulus, Mpsi 

axial loads, kip 

lateral loads, kip 

fiber volume ratio 

in-plane fiber shear modulus, Mpsi 

out-of-plane fiber shear modulus, Mpsi 

matrix shear modulus, Mpsi 

translational spring constant, lb/in. 

torsional spring constant, lb-in./rad 

fiber heat conductivity in longitudinal direction, Btu in./hr in . 2 °F 

fiber heat conductivity in in-plane transverse direction, Btuin./hr in. °F 

fiber heat conductivity in out-of-plane transverse direction, Btu-in./hr in. °F 

matrix heat conductivity, Btu-in./hr in . 2 °F 

torsional moment, kip-ft 

number of fibers per end 

ply delamination strength, ksi 

fiber compressive strength, ksi 

fiber tensile strength, ksi 

matrix compressive strength, ksi 

matrix shear strength, ksi 

matrix tensile strength, ksi 

ply stress 

standard deviation 

ply thickness of skin, in. 

ply thickness of stringer, in. 

void volume ratio 

vector of independent random variables 

fiber thermal expansion coefficient in longitudinal direction, ppm/ °F 

fiber thermal expansion coefficient in transverse direction, ppm/ °F 

matrix thermal expansion coefficient, ppm/ °F 

matrix moisture coefficient, in./in. 

ply misalignment, deg 

in-plane fiber Poisson’s ratio 

out-of-plane fiber Poisson’s ratio 

matrix Poisson’s ratio 

fiber mass density, lb/in. 

matrix mass density, lb/in . 3 
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TABLE I — MATERIAL PROPERTIES AT THE CONSTITUENT 
LEVEL FOR SKIN AND STRINGERS 


Property 

Assumed 

distribution 

type 

Mean* 

Assumed 

uncertainty 

scatter 

E^ Mpsi 

Normal 

3L0 

0.1 

35 

Mpsi 



2.0 



G n2 > Mpsi 



2.0 



G^, Mpsi 



1.0 



v fl2 



0.2 



v f23 



.25 



Of,,, ppm/°F 



-.55 



0 , 22 , ppm/ °F 



5.6 



p f , lb/in. 3 



.063 



N f 

Constant 

10 000 



df, in. 

Normal 

.0003 



Cf, Btu/in. °F 



.17 



Btu* in./hr in. 2 °F 



580 



K q 2 , Btu* in./hr in. 2 °F 



58 



Btu* inVhr in. 2 °F 



58 



S^, ksi 

Wei bull 

400 



ksi 

Wei bull 

400 



E m , Mpsi 

Normal 

.5 



G m , Mpsi 



.185 



v n> 



.35 



«m, ppm/ °F 



42.8 



Pm- lb/in ‘ 



.0443 



C m , Btu/in. °F 



.25 



K„ Btu- inVhr in. 2 °F 



1.25 



S raT , ksi 

Weibull 

15 



S mC’ 

Weibull 

35 



Sms. ksi 

Weibull 

13 



P m ,(iiu r m.)/1% moist 

Normal 

.004 



D m , in. 3 /sec 

Normal 

.002 




“Typical values for graphite-fiber/epoxy-matrix composites at 0.6 fiber 
volume ratio. 


TABLE n.— FABRICATION VARIABLES 
AT PLY LEVEL 


Variable 

Assumed 

distribution 

type 

Mean 

Assumed 

uncertainty 

scatter 

fvr 

Normal 

0.60 

0.05 

VYT 



.02 

.05 

e p , deg 



.00 

.9 (stvd) 

Vsk- m- 



.005 

.05 

w in - 



.02 

.05 
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TABLE HI — UNCERTAINTIES IN STRUCTURAL LEVEL 


Uncertainty 

Assumed 

Mean 

Assumed 


distribution 


uncertainty 


type 


scatter 

It.™, lb/in. 

Normal 

30xl0 6 

0.20 

cTR’ 

K cTO , lb-inirad 
F x , kip 



UxlO 2 

288 

.20 

.05 

F y , kip 



5.76 

.05 

M„, Rip-ft 



576 

.05 
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Probability density function (PFD) 



Figure 2 —Probabilistic design assessment of composite structures. 


Probabilistic design criterion 
Minimum allowable performance X 
Probability (operating performance < X) < 0.001 

Critical performance Y 

Probability (operating performance < Y) = 0.001 


Allowable 
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Free end 



Skin 

cylindrical shell 
[±45/0 2/±45/0 2/±45/0/90] s 

Circumferential frame (horizontal) 
three T-shape frames (10 ft apart) 
[0243 

Stringer (vertical) 

four T-shape stringers (90 deg apart) 

[0 243 

Thickness 

ply thickness 0.005 in. 


Figure 4.— Geometry and composite configuration of stiffened 
composite pipe. 




-O Mxx 








5.2 < F y < 6.4 246 < M xx S 606 


Figure 5. — Probabilistic loading conditions (kip) of stiffened composite pipe. 
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Free end 




Figure 6. — Cumulative distribution function (CDF) and probability density function (PDF) of lateral displacement at free end. 



Figure 7. — Sensitivity factors of lateral displacement at free end at 0.001 failure probability. 
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Sensitivity factors 

Figure 8. — Sensitivity factors for 0.0008 deiamination probability. 
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